Quantum Accelerator Modes from the Farey Tree 
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We show that mode-locking finds a purely quantum non-dissipative counterpart in atom-optical quantum 
accelerator modes. These modes are formed by exposing cold atoms to periodic kicks in the direction of the 
gravitational field. They are anchored to generalized Arnol'd tongues, parameter regions where driven nonlinear 
classical systems exhibit mode-locking. A hierarchy for the rational numbers known as the Farey Tree provides 
an ordering of the Arnol'd tongues and hence of experimentally observed accelerator modes. 
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Precise control of the state and time evolution of quantum 
systems is of critical importance in many areas of physics. 
Tailoring wave packets in Rydberg systems OL producing sin- 
gle photons on demand |2], creating coherent superpositions 
of macroscopic persistent-current states j^t, and controlling 
the production of multiparticle entanglement II, are promi- 
nent examples of "quantum state engineering". Although al- 
most perfect control has been achieved over these systems, 
this can rapidly lose efficiency when influenced by decoher- 
ence or noise. Additionally, generic features of strongly cou- 
pled quantum systems allow for novel and often robust strate- 
gies of quantum control. In such cases, studied in much de- 
tail in the area of quantum chaos, peculiar eigenstates emerge 
which exhibit unexpected localization properties and dynam- 
ics, and are remarkably inert with respect to uncontrolled per- 
turbations. Prominent examples are nondispersive wave pack- 
ets in periodically driven quantum systems J5|], quantum res- 
onances 

m, and stochastic web states iff ritese "strong 
coupling" quantum control schemes rely on underlying classi- 
cal dynamics, which in general is mixed regular-chaotic 
For such a picture to be meaningful it is in general neces- 
sary to approach the semiclassical limit where the classical 
actions accumulated along typical eigenmodes of the system 
are large compared to ft. The quantum system can then "re- 
solve" the intricate phase space structure of classically mixed 
regular-chaotic dynamics, and classical nonlinear stabilization 
phenomena emerge on microscopic scales. 

One of the most ubiquitous of such stabilization phenomena 
in nonlinear classical dynamics is mode-locking. Eigenmodes 
of a periodically driven dissipative system are locked in their 
time-evolution onto the phase of an external drive through a 
nonlinear resonance phenomenon. It occurs in applications 
ranging from frequency-stabilized lasers full to plasma con- 
finement in fusion reactors . An important question is: is 
this necessarily strictly a classical/semiclassical effect? We re- 
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FIG. 1 : Density plot of experimental atomic momentum distributions 
(measured in a frame falling freely with g) after n = 30 pulses as T 
is varied in the vicinity of T% = 133.3yus, from 124.5yus to 142.5/js 
in steps of 0.128/is 1 13]. The labels (a), (b), and (c) indicate (2, 1), 
(3, 1), and (5, 2) quantum accelerator modes, as predicted by the cor- 
responding Arnol'd tongues depicted in Fig. [3] Labels (d) and (e) 
indicate (7, 3) and (8, 3) modes. White lines are the momenta pre- 
dicted by Eq. J3j- Inset shows the evolution of a typical (in this case 
a (1, 0)) QAM as a function of pulse number. The color bar indicates 
the population scale. 



port that features of mode-locking occur at the quantum level 
even far from the semiclassical limit. This was achieved using 
quantum accelerator modes (QAMs) of cold atoms kicked by 
a pulsed standing wave of light orientated along the Earth's 
gravitational field 1 13]. In this Letter we show that the modes 
observed in these experiments can be classified according to a 
number-theoretic construction known as the Farey Tree 1 14]. 

We commence by describing the experimental system 
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which forms the basis of our analysis foil . To create quan- 
tum accelerator modes in the laboratory, laser cooled cesium 
atoms are exposed to a sequence of equally spaced pulses 
from a standing wave of light which is far detuned from the 
nearest atomic transition. Due to the AC Stark effect, the 
atoms experience each ^-function like pulse as a sinusoidal 
potential (spatial period A/2) proportional to the intensity of 
the light. A QAM is characterized by a momentum transfer, 
to a substantial fraction of the atoms, which increases linearly 
with the number of pulses. Figure ^ shows the momentum 
of the atomic ensemble as a function of pulse period T with 
the accelerator modes highlighted by the labelled curves IrHl . 
The inset plots the momentum distribution of a ( 1 , 0) QAM as 
a function of the number of pulses. Note how the width and 
amplitude of the accelerator mode (the peak which moves to 
the left) remain relatively unchanged and the center momen- 
tum increases linearly with pulse number. This can be con- 
trasted with the peak near zero momentum which broadens as 
pulses are applied. Although the stability of the accelerator 
mode from pulse-to-pulse is already somewhat reminiscent of 
mode locking behavior, there are other reasons to think of the 
accelerator mode in this way. 

To see this, we examine the atomic center-of-mass dynam- 
ics using the one-dimensional <S-kicked accelerator Hamilto- 
nian H = p 2 /2m + mgz - h$dW + cos(Gz)] Ym S(t - nT) 1 15]. 
Here z is the vertical position, p the momentum, m the atomic 
mass, g the gravitational acceleration, t the time, G = 4n/A 
llfill . and <pd = Uotp/2h, where Uo is the maximum AC Stark 
shift in the standing wave and t p is the pulse duration. The 
special resonant values of T in the vicinity of which QAMs 
occur experimentally are TV = 2jrfm/{hG 2 ) = ^x66.7jus, with 
t any non-negative integer. Hence, e - 2n((T /TV - 1) is a 
small parameter. Translating to a frame accelerating with g, 
to remove the linear potential, and taking the limit e — > the 
quantum dynamics of the kicked atoms can be modelled by 
the classical map 

J n +\ =J n -Ksm{e n ) - sgn(e)27rO, (la) 
0n+i -On + sgn(e)7„ +1 mod(27r), (lb) 

where sgn(e) is positive/negative if the pulse interval T is 
greater/smaller than TV, and 

9=Gz mod(2w) (2a) 
/„ =I„ + sgn(e)[7rf +j0T- 2nQ.(n + 1/2)], (2b) 

with p/hG = I/\e\ + K = <p d \e\, £2 = gGT 2 /2n, t = 
2jt(T/T(. Note that the deeply quantum mechanical character 
of the atomic dynamics remains hidden in the parametrization 
of J 1 121 through the quasimomentum /?, < (3 < 1 , since the 
limit e — > leading to the classical equations Q leaves the 
finite value of h unaffected. By Bloch theory, subspaces of 
different quasimomenta are decoupled. 

Mode-locking enters the theory of QAMs via Eq. Q which 
also describes the deterministic motion of a periodically 
kicked classical particle on a circle. In this case 9„ and J„ are 
the angle and angular momentum just before the nth kick, K 



is the kicking strength, and O the unperturbed winding num- 
ber. If the classical particle is additionally subject to dissi- 
pative forces, the accessible phase space shrinks, and Eq. Q 
reduces (in the long-time limit) to the sine-circle map 11811 : 
0/7+1 = 0/i - Ksin(6„) - 2nQ., a paradigm in the study of 
mode-locking. If K = and Q is a rational number m/p, 
any trajectory of the sine-circle map returns to its initial value 
(modulo 2n) after p iterations. For < K < 1, mode- 
locking is observed; over a range of Q values around m/p 
(the mode-locking interval) a periodic trajectory with ratio- 
nal winding number m/p persists. This orbit attracts all other 
orbits asymptotically in time, such that finally all have this 
winding number. The widths of the mode-locking intervals 
are exponentially small in p, and increase with increasing K 
up until K — 1 . The regions thus formed in (Q, K) parameter 
space, terminating at K — 0, Q. — m/p, are known as Arnol'd 
tongues 11911 . 

Using this formalism, it is now possible to analyze the dy- 
namics of the QAM in the e — > limit. To begin we look 
for stable periodic orbits in Eq. such that, if {Jo, 9q) is on 
an order p orbit, after p pulses J p mod(27r) = Jq. If the orbit 
is stable, then each of the p points it is composed of is sur- 
rounded by a nonlinear resonance island, set in a chaotic sea, 
where the motion is predominantly regular; the motion in the 
island system approximates that of the periodic orbit. If the 
orbit has winding number m/p, then J v — Jq + 2nm. Thus, 
from Eq. 0, T« (and therefore p„) grows linearly with time. 
The islands travel in momentum, resulting in acceleration. If 
a wavepacket is launched within an island surrounding a sta- 
ble periodic orbit, the acceleration of the corresponding QAM 
obeys 

2n I m\ 

pn^po + n—lQ, «G, (3) 

N \ VI 

precisely as observed in the inset of Fig. [2 Hence, we can 
identify QAMs with nonlinear resonance islands in the clas- 
sical phase space generated by 0, what is just another man- 
ifestation of the general mode-locking phenomenon we are 
describing here. These islands are robust structures, as guar- 
anteed by the Kolmogorov-Arnol'd-Moser (KAM) theorem 
112 OB . A quantum wavepacket initially prepared in the island 
travels with it, and decays only slowly by tunneling into the 
chaotic surroundings. For sufficiently small e, this tunneling 
is exponentially weak, resulting in a stable QAM. More im- 
portantly, due to the KAM theorem the island itself is rather 
inert with respect to perturbations of the Hamiltonian generat- 
ing the map Eq. Q. This robustness is inherited by the QAMs 
and shields them against experimental noise |5]. 

In Fig.|2]we plot a "phase diagram" to represent the regions 
(tongues) where stable periodic orbits with different values of 
(p, rrt) are numerically observed in the (Q, K) parameter plane. 
This plot contains the parameter space explored experimen- 
tally 11311 using values of T in the vicinity of T2 = 133.3/^s. 
Close to K — 0, each of the stable periodic orbit regions is 
wedge-shaped, with its vertex at Q. — m/p, K = 0. Moving to 
higher K inside a tongue, the periodic orbit eventually turns 
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FIG. 2: Arnol'd tongues structure. Different colour identify different 
tongues associated with different quantum accelerator modes. The 
dotted line marks the locus of experimentally explored points when 
T ~ T 2 ; dashed lines bound regions for which (21,8) and (5, 2) stable 
periodic orbits exist, as specified by \K\ > 2n-\fp\Cl - m/p|. 



unstable. A sequence of bifurcations follows, which breaks 
the tongue into fragments. Fragments of different tongues in- 
tertwine and overlap in complicated ways. A tongue may be 
overlapped by others even before breaking, and such overlaps 
persist even at quite small values of K. 

Using a canonical perturbation theory I20H21I E2I1 to deter- 
mine an existence condition for a stable periodic orbit for Eq. 
Q with a given (p, rtt), one obtains, \K\ > 2n^\Q. - rrt/p|. 
In the form of an equality, this equation accurately bounds the 
wedge-shaped (p, m) tongue near its vertex. This is shown for 
the (21,8) and (5,2) periodic orbits by dashed lines in Fig. 
12 Numerical computation and scaling considerations reveal 
the "critical region" where a tongue breaks to be roughly lo- 
cated at kicking strengths K ~ 27rp -3 ^ 2 12211 . Thus, the higher 
the period of an orbit, the narrower the corresponding tongue, 
and the lower the "critical value" of K at which the tongue 
begins to break. 

The parameters corresponding to a specific experiment de- 
termine a point in the phase diagram. If this point is inside a 
tongue then a QAM may be observed. At fixed pulse number 
n, Eq. defines a curve of enhanced population in the (T, p) 
plot (seen in the data of 11311 presented in Fig.0, due to the 
presence of the (p, rtt) QAM. We explore the phase diagram of 
Fig. [2] keeping both <pd and n constant, while varying T. This 
procedure varies e, K and Q. The results of such an exper- 
iment are shown in Fig. [0 The locus of the experimentally 
explored points in the phase diagram is a curve shown by the 
dotted line in Figs.[2]and[3] This curve hits the K = axis at 



Q.' = gGTf/2n, the value of Q. corresponding to the exactly- 
resonant value of the kicking period T = T(. This is 0.3902 
when I = 2. 

Perhaps most remarkably, the values of rrt and p corre- 
sponding to experimentally observable QAM are determined 
by the Farey hierarchy of rational numbers |14]. This rep- 
resentation of rational numbers is a generic feature of mode- 
locking phenomena normally observed in systems with dis- 
sipation. In this hierarchy all rational numbers in [0, 1] are 
constructed as follows: Start from the pair (j, j). At the sec- 
ond level the fraction j - jjj-J- is introduced so that the series 
consists of (j, i, j). On the next level the fractions ~ = y±i 
and I = ^-j-j- are added. This process is continued so that 
if r\ = and — ^ are adjacent irreducible fractions at 
some level, the first rational to be added between them at the 

next level is their Farey mediant r = '"'t'" 2 . At no level can a 

J P1+P2 

rational with a denominator smaller than Pi + P2 be found be- 
tween r\ and r^. At each level the interval [0, 1] is thus divided 
by the Farey fractions into Farey subintervals. As the experi- 
mental line approaches O' in Fig.|3| it successively intersects 
tongues specified by values of (p, m); these values determine 
the observed QAM. The ratios m/p are increasingly close ap- 
proximations to Q'. 

To determine the tongues (and hence QAMs) appearing in 
the experiment, we start from orbits with small p. In Fig.|5K 
the (1,0) and (1,1) tongues are presented (their vertices, at 
K — 0, are outside the boundaries of Fig. 0. These corre- 
spond to the first numbers in the Farey hierarchy. The dot- 
ted line marking the experimental points intersects (within the 
boundaries of the figure) the (1,0) tongue. In the region of 
intersection the stable orbit (1, 0) is found. The corresponding 
QAMs exhibit rapid acceleration, and at n — 30 pulses they 
move beyond the experimental window shown in Fig.^ For 
higher-order QAMs, higher orders of the Farey hierarchy are 
required. At the second level the (2, 1) tongue, shown in Fig. 
|3j}, is introduced, and Q' is in the interval [j, A], The exper- 
imental line intersects the (2, 1) tongue, so a (2, 1) orbit (and 
QAM) is found. Corresponding points are marked by (a) in 
Figs.^andEJ}. The third level (Fig. 03) introduces the (3, 1) 
and (3, 2) tongues. The experimental line intersects both these 
tongues, yet only the (3,1) QAM is observed (intersection re- 
gion marked by (b) in Figs. ^ an d|3p). This is because I is 
further than | from Q.', and so the intersection with the (3, 2) 
tongue takes place in a region where K > 27rp~ 3 / 2 . There are 
only narrow remnants of the tongue, and the corresponding 
stable island is too small for a QAM to be observable. The 
relevant Farey subinterval is now A to j. The construction 
can be continued in similar fashion. In Fig.[3f> the (5, 2) and 
(5, 3) tongues are introduced. Since both lines have a large 
overlap with the (5,2) tongue (regions marked with (c)), the 
corresponding QAM appears on both sides of the resonance. 
Proceeding would ideally produce all the Farey subintervals in 
which Q' belongs. Faint traces of QAMs that lie outside this 
recursion may also be detected, e.g. the white curve (d) in Fig. 
^corresponds to a (7, 3) mode. Note how in the experimental 
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FIG. 3: Operation of the Farey recursion, moving from frames A 
to D, for determining the experimentally observed (2, 1), (3, 1), and 
(5, 2) quantum accelerator modes. The dotted line indicates the locus 
of the experimentally explored points, while the labels (a), (b) and 
(c) mark the regions in which accelerator modes are observed. The 
labels correspond to those used in Fig.Q 

data it is disfavored in comparison with the (8, 3) mode (curve 
(e)), as | is closer than I to Q'. 

This construction demonstrates how the Farey tree classi- 
fies the complex structure of overlapping tongues according 
to those that are most important for the description of QAMs 
observed for a specific value of Q'. Furthermore, as K — > 0, 
the value of j for the QAMs seen in the experiment con- 
verges to Q'. As £1' is determined by the local value of gravity, 
we obtain systematically improving rational approximants of 



g. The underlying classical mode-locking mechanism thus 
renders quantum accelerator modes a robust tool for efficient 
quantum state control, deep in the quantum realm. 
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